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Abstract 

NaFe 2 As 2 is investigated experimentally using powder x-ray diffraction and Raman spectroscopy 
at pressures up to 23 GPa at room temperature and using ab-initio calculations. The results reveal 
a pressure-induced structural modification at 4 GPa from the starting tetragonal to a collapsed 
tetragonal phase. We determined the changes in interatomic distances under pressure that allowed 
us to connect the structural changes and superconductivity. The transition is related to the for¬ 
mation of interlayer As-As bonds at the expense of weakening of Fe-As bonds in agreement with 
recent theoretical predictions. 
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I. INTRODUCTION 


Nature of superconductivity in iron-based superconductors (FeSCs) remain a challenge of 
modern condensed matter physics due to the complicated interplay between their structure, 
magnetism, electronic nematicity, and superconductivity^’^ A generic picture for the early 
discovered sister cuprate system is that superconductivity emerges and develops in a dome 
in the charge carrier concentration versus temperature phase diagram when magnetism of a 
parent compound is suppressed and a narrow spin glass state^ is passed. To study the emer¬ 
gence of superconductivity, one can naturally apply chemical doping, however application 
of pressure has the similar effect in Fe SCs as there are some similarities for the structural 
distortions between pressure and chemical doping.^ In the case of 122 Fe SCs, pressure- 
induced suppression of antiferromagnetic spin fluctuations was proposed to account for the 
disappearance of superconductivity in a so-called collapsed tetragonal (CT) phase.^ How¬ 
ever, the structural modihcation has been suggested to be a key to the understanding of 
the superconductivity of iron pnictides.® A common structural characteristic of all Fe SCs is 
the presence of As-Fe-As layers forming edge-sharing FeAs 4 tetrahedra. It is believed that 
superconductivity emerges from these layers. 

NaFe 2 As 2 crystallizes in the tetragonal (T) ThCr 2 Si 2 -type structure (see 

Fig. 1), belonging to the so-called 122 family of FeSCs . It has been characterized as the 
“missing” member of the 122 system, since initially it was thought that Na, due to its small 
ionic radius, cannot £11 the coordination sphere^’® between the As-Fe-As layers. Moreover, 
NaFe 2 As 2 has a remarkably high superconducting transition temperature Tc of 25 K® in 
comparison with less than 5 K in the case of Cs and Rb within this AFe 2 As 2 family. 

Previously reported structural studies on 122 FeSCs under pressure (e.g. ReF) were 
limited to the cases of A being a divalent alkaline earth or rare earth metal atom, while 
systems with monovalent alkali metal, to the best of our knowledge, have not been studied. 
Although all previous studies suggested that the transformation from the T to the CT phase 
is a universal characteristic of divalent AFe 2 As 2 compounds, a detailed structural study 
on the atomic level and the link between structural modifications and superconductivity 
under pressure is still missing. In this paper we present a combined experimental, using 
x-ray diffraction (XRD) and Raman spectroscopy, and theoretical study of NaFe 2 As 2 under 
pressure. We examine the structural behavior of AFe 2 As 2 superconductors under pressure 
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FIG. 1. (a) ThCr 2 Si 2 -type crystal structure of NaFe 2 As 2 , (b) PbClF-type crystal structure of 
NaFeAs. 

in the case of A being an alkali metal to obtain a deeper insight on the interplay between 
structure and superconductivity under pressure. We present, for the first time, direct ex¬ 
perimental evidence of the formation of As-As interlayer bonds under pressure. According 
to theoretical predictions^^’^^ this is the key parameter determining the correlation between 
structure and superconductivity for 122 SC. 


II. EXPERIMENTAL AND COMPUTATIONAL METHODS 

Fine powder prepared of NaFe 2 As 2 single crystals were loaded into a diamond anvil cell 
(DAC) for the angle dispersive XRD experiments. For the high-pressure Raman measure¬ 
ments, small pieces with a typical dimension of 50 /rm were loaded into a DAC. Neon was 
used as the pressure-transmitting medium for both sets of measurements. Small quantities 
of ruby and gold powder were also loaded, for determination of pressure through ruby lu¬ 
minescence and gold equation of state (EOS), respectively. XRD data were collected at the 
GSECARS (sector 13), using a MAR355 CCD detector. The monochromatic x-ray beam 
(A=0.3344 A) was focused to a nominal diameter of 4 /rm. The images were integrated using 
the FIT2D^^ program to yield intensity versus 26 diagrams. Raman spectra were measured 
in a 135 degree^'^ geometry using the 532 nm line from a solid state laser for excitation. 
Ultralow fluorescence type Ila synthetic diamond anvils were used for Raman experiments. 


3 













Single crystals in the Na-Fe-As system were grown by use of the NaAs flux method. We 
obtained NaAs by reacting the mixture of the elemental Na and As in an evacuated quartz 
tube at 200 °C for 10 h. NaAs and Fe powders were carefully weighed and thoroughly 
ground. The mixtures were put into alumina crucibles and then sealed in iron crucibles 
under 1.5 atm of highly pure argon gas. The crucibles were sealed and heated to 950 °C at 
a rate of 60 °C/h in the tube furnace hlled with the inert atmosphere and kept at 950 °C for 
10 h and then cooled slowly to 600 °C at 3 °C/h to grow single crystals. The shiny crystals 
with typical size of 6x6x0.2 mm can be easily cleaved from the melt. 

Our theoretical calculations were performed using density functional theory within the 
generalized gradient approximation^® and using the PAW method^® as implemented in the 
VASP code.^^ The PAW potentials used had [Be] core (outermost radius 2.2 a.u.) for Na 
atoms, [Ar] core (radius 2.3 a.u.) for Fe atoms, and [Ni] core (radius 2.1 a.u.) for As atoms. 
We used the plane wave kinetic energy cutoff of 360 eV and 6x6x3 F-centered meshes for 
Brillouin zone sampling. Electronic optimization was done self-consistently with a threshold 
of 10“^ eV/cell, and structure relaxation proceeded until changes in the enthalpy were below 
10“^ eV/cell. A ferromagnetic conhguration with all Fe atoms given a starting magnetic 
moment of 2 Bohr magnetons was used. The starting magnetic moment plays little role (only 
the magnetic symmetry is important) and was optimized during electronic and structural 
relaxation. 


III. RESULTS 

X-ray diffraction results and analysis 

Figure 2 shows selected XRD patterns of NaFe 2 As 2 obtained at various pressures up to 
23 GPa. All the observed peaks in this pressure range can be properly indexed according to 
the tetragonal ThCr 2 Si 2 -type structure. On the other hand, an unusual pressure behavior 
of the Bragg peaks related with the a-axis can be clearly seen at low pressures. This is 
highlighted by the merge of the (114) and (200) peaks at 26 ~ 10° and the splitting of the 
peak at 26 ^ 7° into its components, (110) and (103). We observe a decrease of the 26 
position of the (200) peak and a practically constant position of the (101) peak, while the 
26 positions of all other peaks [e.g. (002)] shift towards higher values. 

In order to determine the positional parameters and the interatomic distances we have 
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FIG. 2. X-ray diffraction patterns of NaFe 2 As 2 at various pressures. The asterisks denote peaks 
from Ne. 

performed a detailed Rietveld refinement of the XRD patterns using the GSAS program.^® 
Typical rehned prohles are shown in Fig. 3 at 1.0 and 5.1 GPa. The only free positional 
parameter of this structure is the 2 : position of As atoms, Wyckoff position (WP) 4e (0, 0, 
z), since positions of Na and Fe atoms are hxed. Our two-dimensional (2D) images reveal 
an almost perfect ring shape of reflections, mainly due to the very fine powder used in this 
study. In addition wide angle dispersive XRD measurements helped us to overcome any 
intensity redistribution caused by preferred orientation. Moreover, XRD patterns have been 
collected from several spots of the sample at the same pressure, and the results were aver¬ 
aged. The agreement between experimental and theoretical interatomic distances provides 
an additional confidence of these results. 

From XRD data of NaFe 2 As 2 we have obtained the lattice parameters and the unit cell 
volume as a function of pressure and the results are shown in the plots of Fig. 4. The pressure 
evolution of the lattice parameters and the c/a ratio changes at a critical pressure Pc ~ 4 
GPa. Below P^, a-axis increases while both the c-axis and the c/a ratio decrease rapidly. 
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FIG. 3. Rietveld refinements for NaFe 2 As 2 at (a) 1.0 GPa and (b) 5.1 GPa. Respective 2D x-ray 
diffraction images are shown in insets. 

In contrast above a-axis starts to decrease and c-axis becomes much less compressible. 
Consequently the c/a ratio shows a gradual decrease. The different behavior of c-axis, below 
and above Pc-, can be viewed as a modihcation of the initial tetragonal phase (T) to the so- 
called collapsed tetragonal (CT) phase [e.g. Ref.^®). We have htted the pressure-volume 
data by the third-order Birch equation of state^° and determined the bulk modulus Bq and 
its hrst pressure derivative Bq (Fig. 4). The CT phase is much less compressible compared 
to the T phase as Bq increases by a factor of 5. No apparent volume discontinuity is observed 
at this isostructural phase transition. 

The results from the theoretical predictions are also plotted in Fig. 4 for both a magnetic 
(M) and a non-magnetic phase (NM). In the case of the M phase, a ferromagnetic conhgu- 
ration has been used and the energy difference from the antiferromagnetic conhguration is 
very small. It is well known that density-functional calculations have a systematic error of 
several GPa in the equation of state, and a simple constant shift usually brings computed 
equations of state into close agreement with experiment.Here, the theoretical pressures 
were corrected by 4 GPa. With this correction, an almost perfect agreement between ex- 
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FIG. 4. Pressure dependence of (a) lattice parameters and (b) cell volume of NaFe 2 As 2 . The inset 
shows the pressure dependence of the c/a ratio. Experimental data are shown with solid symbols 
and theoretical predictions (shifted by 4 GPa to higher P, see text) with blue (M-phase) and red 
(NM-phase)lines, respectively. The solid black curves in (b) are the third-order Birch-Murnagan 
EOS fits for the T {Bq=22 GPa, B'= 4) and GT (5o=110 GPa, B'= 4.5) phases. 

perimental and theoretical values (for both phases) can be clearly seen from Fig. 4. It 
is noteworthy that the M phase, although it has higher volume at low pressures, is more 
stable than the NM one, due to the magnetic interactions which result in a (slightly) lower 
enthalpy up to well above Pc- Moreover, both phases exhibit a T to CT transformation at 
the same pressure. In the magnetic case (which corresponds to observations of Ref.^) our 
calculations suggest that CT phase remains ferromagnetic at T = 0 K. Since theory predicts 
the same phase transformation in magnetic and nonmagnetic cases, we can conclude that 
the structural behavior is not driven by magnetic interactions. 

Raman spectroscopy 

Raman spectra of NaFe 2 As 2 under pressure show three out of four Raman-active zone- 
center modes predicted from group theory:^^ Aig (As(z)), Rig(Fe(z))and i?g(As-Fe(xy)and 
Fe-As(xy)) (Fig. 5(a)). Aig and Big modes correspond to the displacement of As and Fe 
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atoms along the c-axis, respectively. The two Eg modes correspond to mixed lateral As 
and Fe displacements inside the As and Fe layers. The ambient pressure frequencies are in 
very good agreement with previous studies^^ with the exception of the very low intensity 
high frequency F^g(Fe-As(xy)) mode, which was not observed. All Raman modes show 
normal mode behavior he. increase of the frequency with increasing pressure in phase T 
(Fig. 5(b)). At 3 GPa, a critical pressure which is in agreement with that determined from 
XRD measurements, there is an apparent change in the pressure slope of Raman modes 
(see Table I), accompanied by an increase of the linewidth. Remarkably, the Big{Fe) mode 
shows a softening in CT phase. The mode Griineisen parameters {'Jt) determined using the 
experimental results of this work are shown Table I. The jt parameters of modes of the 
T phase are common for materials with mixed covalent-ionic bonding (such as within the 
As-Fe-As layers), due to the presence of weaker interlayer bonds, which experience a larger 
compression. In the GT phase, the Aig (As) mode exhibit a higher 'Jt which almost doubles. 
These changes suggests an alteration of this mode character from intralayer to interlayer , 
which results in an increase of the contraction of the respective bond. Similarly, we hnd the 
doubling of the yr of the intralayer Eg (As-Fe) mode. It is the most intriguing that the yr 
of the Big (Fe) mode becomes negative, in GT phase suggesting weakening of bonds under 
pressure. 
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FIG. 5. (a) Raman spectra of NaFe 2 As 2 at various pressures, (b) Raman peaks frequencies plotted 
against pressure. In (b) lines are the results of linear fits. 
















TABLE I. Assignments, ambient conditions frequencies wq (in cm“^), slopes {dio/dP)T (in 
cm“^GPa“^) and Griineisen parameters /d{inV)\T of the Raman modes of the T and 

GT phases. 



T 

GT 

Mode 

CJO 

{d(jj/dP)T 

It 

UJQ 

{d(jj/dP)T 

It 

E9 

110.7 

5.3 

1.05 

119.8 

2.4 

2.2 

Alg 

180.7 

8.6 

1.05 

197.8 

3.4 

1.9 

Big 

200.9 

17.0 

1.86 

259.2 

-1.7 

-0.72 


IV. DISCUSSION 

Our x-ray diffraction results reveal the behavior of the As-As interlayer distance in com¬ 
parison with the Fe-As intratetrahedral distance (see Fig. 6). Figure 7 shows the experi¬ 
mentally and theoretically determined values of these distances under pressure. An almost 
perfect agreement between the two sets of values is observed. The As-As interlayer distance 
decreases rapidly with pressure up to 4 GPa. This decrease is almost half of that of c-axis 
and it accounts for the whole c-axis contraction while the thickness of As-layers (vertical 
distance between As-As layers of the same As-Fe-As block) remains almost constant or even 
increases. Shortening of c-axis is driven exclusively by the distance between the tetrahe- 
dra blocks {c=2dAs-As(c)+‘^dAs-Fe-As(c))- The interlayer As-As and the intralayer Fe-As 
distances become practically equal above 4 GPa, i.e. after the transformation to the GT 
phase. The As-As distance above 4 GPa takes a value (2.4 A in experiment and 2.6 A in 
theory) which is consistent with the formation of As-As bond (covalent radius of As is 1.2 
A23). Thus, we conclude that a major change in bonding between As atoms occurs under 
pressure. It can be viewed as a change from 4-fold to 5-fold coordination with formation of 
interlayer As-As bonds (see Fig. 6(b)). This conclusion is further justihed by the doubling 
of the 7 r of the Aig mode after the phase transition. In contrast the decrease of the of 
the Big mode indicates weakening of the respective Fe bond. 

Previous experimental studies on AFe 2 As 2 (A=Eu,^° Ba,^^ Sr,^® Ga^®) compounds, reveal 
a common structural trend under pressure. This is an increase of a while both the c and 
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FIG. 6. (a) ThCr 2 Si 2 -type crystal structure of NaFe 2 As 2 showing the different interatomic distances 
and (b) ThCr 2 Si 2 -type crystal structure of NaFe 2 As 2 with As atoms in 5-fold coordination. 

c/a ratio decrease rapidly below Pc, and a decrease of a and mnch less compressible c 
above Pc- In the light of this observation, Uhoya et snggested that the transformation 
from the T to the CT phase is a nniversal characteristic of AFe 2 As 2 componnds, while 
the valne of Pc depends on the divalent metal. Onr stndy reveal a similar behavior for 
single valence metal snggesting that this nniversality is also independent of the valence of A 
metal. Transformation to the CT phase with farther increasing pressnre always snppresses^^ 
or diminishes snpercondnctivity.^® In a stndy of CaFe 2 As 2 nnder pressnre^® it has been 
snggested that one of the main strnctnral differences between T and CT phases is the valne of 
the As-Fe-As tetrahedral angles. Both tetrahedral angles change abrnptly dnring the phase 
transition from a valne very close (T) to a valne far away (CT) from the ideal tetrahedral 
angle (109.5°). Since decreases with increasing deviation from the ideal tetrahedral angles, 
it has been conclnded that this is crncial for snpercondnctivity althongh the direct link and 
the nnderling mechanism which correlates tetrahedral angles and snpercondnctivity is not 
clear. 

Recent theoretical stndies^^’^^’^® were focnsed on the change in interatomic distances 
nnder pressnre, mainly on the interlayer As-As distance (noted as As-As (c) in Fig. 6). This 
distance has been predicted to decrease abrnptly with pressnre, reaching a valne that is 
very close to the As-As covalent bond distance in the CT phase. Conseqnently it has been 
proposed that the transition to the CT phase is indnced by formation of a direct As-As 
interlayer bond and a weakening of the in-plane Fe-Fe^^ and Fe-As bonds.It is noteworthy 
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Pressure (GPa) 

FIG. 7. Pressure dependence of selected interatomic distances of NaFe 2 As 2 : As-As interlayer along 
the c axis and Fe-As inside the FeAs 4 tetrahedra. Experimental and theoretical values are shown 
with black symbols and red lines respectively, the dashed lines are guides for the eye. The inset 
shows the pressure dependence of the two distinct tetrahedral As-Fe-As bond angles. The ideal 
tetrahedral angle is indicated by the horizontal dashed line. 


that in the case of 111 superconductors (e.g. Ref^^) the Fe-As-Fe layers are separated by 
two layers of cations (see Fig. 1). Thus, an interaction between As atoms of different layers 
is hindered and only intralayer modifications can affect superconductivity. 

The spin-state of Fe is the key parameter which controls As-As bonding and consequently 
the lattice parameters. Under compression the magnetic moment of Fe decreases and con¬ 
sequently, the strength of the Fe-As bonding (see Fig. 4 of Ref.^^). This is also evident from 
our calculations, see Fig. 8. The results of our calculations suggest that the CT phase is 
magnetic i.e. the Fe magnetic moment is not vanished but suppressed after the transforma¬ 
tion decreasing the strength of the Fe-As bonding in agreement with softening of the Big 
Raman mode (Fig. 5). The structural changes determined in this paper (Fig. 7) clearly 
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support this scenario. A direct consequence of the As 5-fold coordination is that the sys¬ 
tem loses its two-dimensionality (Fig. 6(b)), and with it, superconductivity. Moreover the 
weakening of the in-plane Fe-Fe bonds^^ in the CT phase also decreases superconductivity. 



FIG. 8. Calculated magnetization per unit cell of NaFe 2 As 2 as a function of pressure. 

V. SUMMARY 

From our data we hnd that the tetrahedral angles change continuously with pressure 
(inset Fig. 7) in a manner very similar to that determined previously in CaFe 2 As 2 ®’^® and 
also predicted by theory. Dnring the initial compression the a-axis expands and the c-axis 
rapidly decreases dne to the development of the As-As interlayer bonding and this keeps 
Fe-As distances almost constant. In addition the tetrahedral angles are close to, or even 
approach, the ideal valne. Above Pc the tetrahedral angles take valne away from the ideal 
because a-axis stops to increase to mediate the c-axis decrease. So, the correlation between 
the deviation of tetrahedral angles and Tc turns ont to be a “side effect” of the change of As 
atom position with increasing pressure which primary affect . In conclusion, our study 
reveals, for the first time, the strong experimental evidence snpporting the theoretically sng- 
gested mechanism controlling the structnral behavior of AFe 2 As 2 componnds nnder pressnre. 
Moreover the concomitant experimental and theoretical study of NaFe 2 As 2 provides, miss¬ 
ing up to now, a direct link between strnctnral characteristics and snpercondnctivity under 
pressure. 
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